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Abstract 

We address the "inverse problem" for discrete geometry, which consists in determining whether, 
Q ' given a discrete structure of a type that does not in general imply geometrical information or even 

a topology, one can associate with it a unique manifold in an appropriate sense, and constructing 
the manifold when it exists. This problem arises in a variety of approaches to quantum gravity 
that assume a discrete structure at the fundamental level; the present work is motivated by the 
semiclassical sector of loop quantum gravity, so we will take the discrete structure to be a graph 
and the manifold to be a spatial slice in spacetime. We identify a class of graphs, those whose 
vertices have a fixed valence, for which such a construction can be specified. We define a procedure 
designed to produce a cell complex from a graph and show that, for graphs with which it can be 
carried out to completion, the resulting cell complex is in fact a PL-manifold. Graphs of our class 
for which the procedure cannot be completed either do not arise as edge graphs of manifold cell 
Q\ decompositions, or can be seen as cell decompositions of manifolds with structure at small scales 

(in terms of the cell spacing). We also comment briefly on how one can extend our procedure to 
K*" more general graphs. 



cr 

bj:) 



PACS numbers: 04.60.Pp, 02.40. Sf, 05.90.+m. 



*Electronic address: luca@phy.olemiss.edu 
''^Electronic address: 'corichi(9)matmor . unam. mx 
^Electronic address: .owinklerQperimeterinstitute.ca 



1 



I. INTRODUCTION 



One of the most important trends in the development of our current understanding of space- 
time and gravity has been the decrease in the number of background, non-dynamical struc- 
tures used in formulating the theory. For example, in the canonical approach to quantum 
gravity [H El |3] one starts just with a background different iable manifold, interpreted as 
space, and builds a diffeomorphism-invariant theory without any additional structures, such 
as a preferred metric or coordinate system, on it. A further step in the direction of getting rid 
of "ideal elements" would be to formulate the theory without using a manifold as part of the 
background structure . In addition to a greater conceptual simplicity, we would then have 
a more flexible theory in which even topological properties of spacetime could be seen as dy- 
namically determined. The replacement of the notion of manifold points by a different basic 
structure should set in at a characteristic length scale representing the fuzziness of quantum 
spacetime, and help eliminate the divergence problems that appear in classical spacetime. 
This viewpoint has often been advocated over the past several decades (for example, it can 
be found in an implicit manner in reference [5]), and motivates our work. 

In loop quantum gravity [H El E] the basic building blocks for quantum states are graphs 
embedded in a given spatial manifold; this background manifold is required for the definition 
of the full (kinematical) Hilbert space and of some of the operators. One can envision however 
a version in which abstract graphs and states on them are the fundamental objects, while 
the manifold is an emergent concept; indeed, a variant of loop quantum gravity of this type, 
algebraic quantum gravity, has recently been proposed [B]. In such a formulation canonical 
quantum gravity, until now a theory with strong indications of discreteness, would join a 
number of other theories that are fundamentally discrete, but with the advantage that its 
structure and relationship with classical general relativity are considerably better understood 
than those of most other theories. 

If the manifold is not present from the beginning, however, the question of deciding 
whether a given quantum gravity state is semiclassical, i.e., whether it approximately de- 
scribes a classical geometry, acquires a new aspect with respect to other approaches. The 
issue now is not just whether observables defined on the underlying space are peaked around 
values of the corresponding classical quantities, but whether the underlying structure itself 
resembles a classical space. In some other theories, the way in which a manifold can be 
associated with the discrete structure is in principle straightforward. For example, in the 
case of Regge calculus or dynamical triangulations [71 IH], the discrete structure is a simpli- 
cial complex, which is already a topological space; although one may decide not to use this 
information, and define for example an effective dimensionality using other arguments, each 
fundamental configuration is directly a PL-manifold. The discrete structure of loop quantum 
gravity (as well as that of other approaches such as causal sets [9J) does not have this feature. 
Some simple graphs naturally suggest manifolds associated with them (for example, a cubic 
graph and IR*^), but using only those graphs is overly restrictive (at least if all edges of the 
graph are "active" in the quantum state; and if they are not, the graph becomes effectively 
different). It is true that in three or more spatial dimensions one can always embed any 
graph in a manifold; in fact, in any manifold. However, in order for the theory to really 
be background-free, the manifold, if any, must be determined by intrinsic properties of the 
graph itself; at scales larger than that of the graph, it must not contain additional structure 
not "sampled" by the graph, and at smaller scales it should have no structure at all. 
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Therefore, the question can be loosely phrased as follows: Is there a procedure by which, 
given a graph, one can determine whether there is a manifold that approximates it in the 
above sense, and possibly construct such a manifold if one exists? This problem has been 
called the inverse problem for discrete geometry [S], and its resolution is important for the 
semiclassical sector of the theory in its manifold-free version. The corresponding problem 
in the causal set approach has been addressed in work on causal set kinematics, and some 
limited results have been obtained [TT]. but we are not aware of results in loop quantum 
gravity applicable to the problem as formulated above. However, some work on graphs has 
been done motivated by other theories, aimed at determining the effective dimensionality of 
a graph from scaling considerations between volumes and lengths, as estimated in the graph, 
and from the dimensionality of simplices that can be defined from graph edges [12]. 

Let us comment on the extent to which such proposals will be incorporated in our work. 
The fractal-like definition of dimension based on scaling is a very general one, with the 
advantage that it can always be applied, and the results interpreted in terms of whether 
a graph can be meaningfully assigned a dimension or not. However, it has the limitation 
that it is a statistical definition. In itself, this is reasonable, since one assumes that the 
graph vertices are to be interpreted as uniformly distributed points in the manifold, but it 
also means that the proposal does not use the structure of the graph in detail, nor does it 
provide more information on the manifold in addition to dimensionality, and it is difficult to 
imagine an extension that would. Our point of view is that, ultimately, the most useful way 
of characterizing the dimensionality and topology of a discrete structure will be determined 
by physical arguments, based on the dynamics of such structures and the effective matter 
fields we observe as coupled to them. Our current understanding of the theory does not 
allow us to do this yet, and for the time being, our goal will be to propose an approach that 
uses the full structure of the graph. However, because we will also think of graph vertices 
as uniformly distributed, our proposal will in effect incorporate the results a fractal-like 
definition would give, in situations where both are applicable. 

The definition based on simplices does not have the limitation just mentioned, and it would 
appear to be a good one for our purposes. In fact, a procedure obtained as an extension of 
the simplex-type definition of dimension could lead to good approximating manifolds, if used 
with the right type of graphs. We will not follow this approach just because the physical 
interpretation we follow for graphs in loop quantum gravity suggests that we identify them 
with sets of edges of a different type of cell complex; simplicial complexes can also be obtained 
in our case, but only as a result of a further construction rather than directly. 

This leads to an important point about the inverse problem, which is that we do not 
expect to find a procedure which always produces a "correct" answer, independently of the 
type of graph it is applied to and of the role the graph plays in the theory. In loop quantum 
gravity in general, graphs are not assumed to have any special topological properties, and 
most abstract graphs don't need to "look like a manifold". However, the type of metric 
information a graph can carry in quantum geometry does depend on its intrinsic properties, 
so some graphs are more interesting than others for the semiclassical sector, and in our 
approach to the inverse problem we will limit ourselves to a class of graphs that is well- 
motivated (and large enough) in this context. In sections [H] and HI we spell out what we 
mean by a manifold that approximates a graph, motivate our choice for the class of graphs we 
use and discuss their properties, and write down a precise statement of the inverse problem 
for our class of graphs. In section [TV] we use those properties to address the inverse problem. 
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II. STATEMENT OF THE PROBLEM 



In this section we translate the vague wording of the inverse problem used in the Introduction 
into a precise statement, and specify the class of graphs we will focus on; we begin with one 
definition of what it means for a manifold to approximate a graph 7. The general idea is to 
look for a manifold M in which 7 can be embedded so that there is a "good match" between 
the image of 7 and M, using a cell complex as an intermediary: 

Manifold tiling graph: We will say that a graph 7 is a tiling graph for a PL-manifold 
M if there exists a cell complex VL which is PL-isomorphic to M and such that the 
set of 0-dimensional and 1-dimensional cells (the "1-skeleton" ) of is the graph 7. 

Notice that our definition does not specify the type of cell complex one obtains, for example 
whether it is a triangulation or not; this will depend at least in part on the type of graph we 
consider. Also, it is in principle possible that more than one, inequivalent manifolds can be 
found that fit the definition with the same graph. Using this definition, the inverse problem 
for graphs acquires the meaning of determining whether a given graph is manifold-tiling; but 
before we give a procedure for addressing this question, we need to specify the set of graphs 
we are interested in, discuss the possible manifold non-uniqueness, and comment on whether 
the definition captures what we want from a physical point of view. 

To identify a suitable class of graphs to work with, we appeal to basic facts in geome- 
try and loop quantum gravity. Having chosen cellular decompositions to mediate between 
graphs and manifolds, we notice that in a generic geometry without symmetries, the only 
cell decompositions that can be covariantly defined are the ones that result from a random 
process on the manifold. The most natural ones use uniformly random (Poisson) point pro- 
cesses, in which each point becomes either a vertex (in a Delaunay triangulation) or a cell 
(in a Voronoi complex). In loop quantum gravity, graph edges are the elements on which the 
holonomy variables are defined, with curvature associated with closed loops, while quantum 
geometry results indicate that vertices of the graph on which holonomies are defined are the 
basic elements of spatial volume [13], and that a necessary condition for them to give a non- 
vanishing contribution in a Z)-dimensional manifold is that they be at least [D + l)-valent. 
Now, in combinatorial geometry, extensive geometric information is usually associated with 
elements of the Voronoi complexes; in particular, volume is associated with Voronoi vertices, 
which generically are precisely {D + l)-valent. On the other hand, scalar curvature, for ex- 
ample, is associated with codimension-2 simplices (a fact commonly used in Regge calculus 
and dynamical triangulations), or Voronoi 2-cells. This agrees with common usage in gauge 
theories, where holonomies and connections correspond to curvature and are also associated 
with (fluxes through) Voronoi 2-cells [HI [15] . We conclude that the most meaningful graphs 
for semiclassical loop quantum gravity are the edge graphs of Voronoi cell complexes. 



As we will see in section III, this implies that, when constructing a D-dimensional mani- 
fold, we will only use {D + l)-valent graphs, which may at first sound like an unreasonable 
restriction. If our goal was simply that of modeling classical geometries using discrete graphs, 
then this choice would not be any more restrictive than, say, using simplices and triangula- 
tions, since in the latter all Z)-cells have exactly D + 1 vertices and facets. Voronoi complexes 
and {D + l)-valent graphs are not the only way to discretize manifolds, but they are the 
simplest discrete structures among which we can prescribe a covariant way of picking a ran- 
dom element as a discretization for a given Riemannian manifold. Other graphs that are 
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often used, such as those defined by regular lattices, are not only much more special than 
the ones we consider, but their choice also requires a high degree of symmetry on the mani- 
fold, since we know of no covariant procedure, random or not, that will produce them in a 
manifold without symmetries. Furthermore, even when they can be defined, the long-range 
order of regular lattices introduces additional large-scale structure on the manifold. Thus, 
in principle, to justify the use of such graphs to discretize a manifold with symmetries, one 
would have to argue either for something like a symmetry breaking mechanism, or that the 
extra structure does not alter the physical conclusions one is interested in. 

On the other hand, from the point of view of a theory with basic variables that are 
more general than manifolds, our choice of graphs does amount to a real restriction. We 
will consider this choice as akin to using coherent states in ordinary quantum theory; the 
latter are not the most general semiclassical states, but they are a good model for many 
aspects of semiclassical physics and a good starting point for the analysis of more general 
ones. In our case, physical loop quantum gravity states may have to be based on more than 
a single graph, and quantum fiuctuations in the semiclassical sector may lead us to include 
graphs with obstructions to embeddability in a manifold at some scale, since most graphs 
are not manifold-tiling. Such obstructions may show up in different ways; in cases where the 
obstructions can be attributed to just a few cells, one possible distinction is between local 
ones, in which vertices of would-be neighboring cells are not connected in the appropriate 
way, and non-local ones, in which some graph edges connect what would otherwise be distant 
cells in the complex. Non-locality has long been considered an inescapable aspect of quantum 
gravity, and it has recently been pointed out [12] that it is likely to the linked to the presence 
of non-manifold-like aspects in the graphs used in loop quantum gravity; simulations with a 
toy model for the dynamics of graphs [17] confirm this expectation. We will start commenting 
in section |V] on situations in which the relationship between 7 and Q can be relaxed to some 
extent, to deal with certain types of obstructions to strict embeddability. 

Our goal for the rest of this paper is to find a procedure by which, given a. {D + l)-valent 
graph, one can determine whether it is a tiling graph for a D-dimensional PL-manifold M. 
We will see that there is a simple constructive procedure by which the manifold can be 
found whenever it exists and does not have structure on (combinatorial) scales of the order 
of, or smaller than, the cell size of the complex; by construction it will then be unique. The 
argument is based only on concepts from combinatorial topology, as opposed to combinatorial 
or differential geometry, but one could try to go beyond those results and ask whether 
the PL-manifold has a differentiable structure, if so whether it is unique, and there is a 
Riemannian metric gab such that the graph is the set of edges of a Voronoi complex based 
on a set of points {pi} in M. It turns out that, at least in 2 and 3 dimensions, there is a 
unique differentiable structure for a tiling graph of fixed valence, and suitable metrics can 
be found. In this connection notice that the graph-topological manifold correspondence 
is many-to-one; given any manifold, we can obtain many "statistically equivalent" {D + 1)- 
valent tilings of it simply by choosing different sets of random points, and many "statistically 
inequivalent" ones by using different metrics. Thus, one may be able to use at least part of 
the additional information carried by a specific graph to determine a class of metrics that 
could have produced the corresponding cell complex. We will sketch one way of obtaining one 
representative of such a class, but if one is interested only in quantities such as curvature 
averaged over microscopic regions, statistical techniques are available that do not require 
knowledge of the exact metric, as explained in more detail elsewhere [18] . 
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Finally, let us comment on the mathematics literature regarding the relationship between 
graphs and cell complexes, and the extent to which it applies to our problem. There is a set 
of results on polytopes in combinatorial topology [19j that is closely related to our subject. 
A polytope is the convex hull of a finite set of points in H" (an n-dimensional version of 
a polyhedron). A cell complex in general is not a polytope, but a cell complex that is 
homeomorphic to a sphere is the boundary complex of a polytope in one higher dimension. 
Thus, many results about ra-polytopes become results about cell complexes homeomorphic 
to S"~^, and some results about local properties of polytope faces can be translated into 
results about all {n — l)-dimensional cell complexes (think of n — 1 as D). 

First of all, a graph can be the edge graph of a (D + l)-polytope (or of its boundary 
Z)-complex) only if it is [D + l)-connected,^ which implies that its vertices are at least 
{D + l)-valent. It is also known that a graph can be realized as the edge graph of a 3- 
polytope boundary complex iff it is planar and 3-connected; and that, if a polytope is simple 
(i.e., each vertex is on the boundary of exactly n oi D + 1 cells of codimension 1), then its 
edge graph determines uniquely the polytope. As a whole, there is some overlap between 
the mathematical results we are aware of and ours, but they also complement each other. 

III. VORONOI CELL INCIDENCE RELATIONS 

The purpose of this section is to obtain a set of relations among cells of different dimen- 
sionalities, that are satisfied by generic Voronoi complexes in any manifold; the arguments 
leading to those relations use a Riemannian (i.e., positive-definite) metric on the manifold, 
but the relations themselves are independent of what metric was used. We will also take 
the opportunity to point out with a couple of examples what might happen, in terms of the 
Voronoi complex, when the Riemannian manifold has structure on small scales compared to 
those at which it is being sampled by the cells of the complex. 

Consider a locally finite set of points {pi} at arbitrary locations in a D-dimensional 
manifold M with a Riemannian metric gab-, where the local finiteness condition means that 
every point p & M has an open neighborhood containing only a finite number of p^'s. Of the 
different procedures known that produce graphs in Riemannian geometries, the one based on 
the Voronoi construction (together with its dual Delaunay triangulation) is the most natural 
and, as mentioned in section [ll| the choice is motivated both by simplicity (it gives the 
simplest graphs on which loop quantum gravity states can be defined that encode volume 
information on {M,gab) — see, e.g., references [121 [IHl ED]) and by analogy with the use of 
graphs in combinatorial geometry and discretized gauge theory. 

In the Voronoi construction, the points {pi} act as "seeds" for a partition of M into 
regions Ui, one for each pi (see figure [T]); the Ui and their boundaries then define a cell 
complex Q that is homeomorphic to M. The cell Ui is defined as the set of all manifold 
points which are closer to pi than to any other seed, with respect to the metric gab, while 
the separating {D — l)-dimensional face Uij between Ui and uj is made of manifold points 
that are equidistant from pi and pj, and sets of points that are equidistant from more than 2 
seeds define cells of lower dimensionality, if the seeds are at generic locations. In particular. 



A graph with fc + 1 or more vertices is fc-connected if the subgraph obtained by removing any set of A: — 1 
vertices, together with the edges that end at those vertices, is still connected. 
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FIG. 1: Example of tiling of 2D flat Euclidean space with a random Voronoi complex. The dots are 
the randomly located "seeds" that were used to generate the cell complex. Some vertices appear 
to be 4-valent because of the limited resolution of the figure. 

the vertices of the complex are the manifold points that are equidistant from D + 1 seeds. 
Once we have obtained a Voronoi complex decomposition of M, the last step is trivial: to 
obtain a graph (the "Voronoi graph"), just retain the vertices and edges of the cell complex 
and discard all higher-dimensional cells. 

If the Voronoi complex Q was obtained from points at generic locations in the manifold, 
then the following are satisfied. 

Voronoi cell incidence properties: For any Z-dimensional cell a; in a random Voronoi 
complex Q, the number of fc-cells (with I <k < D) that have uj on their boundary is 

(3.1) 

Thus, each vertex is generically shared by A'^iio = D + 1 edges (for example, in two dimensions 
all vertices are trivalent, while in three dimensions they are all four-valent), and by A^dio = 
D + 1 cells of dimension D, each of which is identified uniquely by specifying which of the 
D + 1 edges at that vertex is not on its boundary. Similarly, each codimension-n cell is on 
the boundary oin + 1 cells of dimensionality one unit higher. To prove these relations, recall 
that a cell of dimension / is the locus of points that are equidistant from a certain number 
n of seeds; since this is equivalent to imposing n — 1 equations on those points, generically 
the dimensionality of the cell is / = D — (n — 1), and therefore n = D + 1 — I. If is on 
the boundary of a cell of higher dimensionality /c, that cell is the locus of points that are 
equidistant from some number m < n of seeds chosen among the n that define u, satisfying 
m = D + 1 — k. The number of such higher-dimensional cells is the number of ways in which 
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m points can be chosen among the n, i.e., 





(3.2) 



Notice that the statement that the relations (3.1) hold for points at generic locations means 



that sets of locations for which they do not hold constitute a subset of measure zero of all 
possible sets of locations. If we used the volume element of the metric on M to sprinkle 
the points uniformly at random, then those relations would hold with probability 1 for the 
resulting set of points, but we will keep our results more general by not invoking the notion 



of random sprinkling here. The relations (3.1) will be very important below, in particular 



because, as we will see, any cell complex f2 satisfying them is equivalent to a PL-manifold. 
They are, of course, not necessary for a cell complex to be a PL-manifold (for example, any 
barycentric decomposition of a cell complex that does satisfy them will produce one that does 
not), but the fact that a cell complex Q does satisfy them provides additional information 
that may be used towards showing that Q can be interpreted as the result of the Voronoi 
construction from some metric. 

On the other hand, for each cell u, the number of lower-dimensional cells that are on its 
boundary is not a fixed number but a set of variables, which tells us how many neighbors 
uj has. The specific set of values that these variables have is what distinguishes one Voronoi 
complex from another and can be used, using appropriate statistical arguments, to estimate 
properties of a metric that is compatible with a given Voronoi complex. 

Let us comment briefly on what might happen if the manifold has structure on small 
scales. Intuitively, we would say that (M, gab) has no structure below the length scale i if 
neither the topology (e.g., from the sizes of non-contractible homotopy generators) nor the 
curvature of the manifold (from the values of curvature invariants) can be used to determine 
lengths of the order of i or smaller. In a Riemannian manifold, we can capture this idea 
by saying that, for every p G M, the ball B£{p) of radius i around p is a convex normal 
neighborhood. Then, if we sprinkle points with density p in a region of a manifold M which 
does have structure below the length scale i = p~^^^, the "cell" Uj around a seedjOj in 
that region may not be homeomorphic to a ball in IR^. In the example of figure |2| the 
Voronoi complex is not a cell complex, strictly speaking, and the graph one obtains from it 
is disconnected (but it is homeomorphic to M, while the triangulation dual to the Voronoi 
complex is 1-dimensional in this region of the manifold, and does not tile it — this is the sense 
in which our Voronoi complexes are not always equivalent to triangulations, as mentioned 
in the Introduction). If the seed spacing was somewhat smaller, but still comparable to 
some characteristic length of the manifold, such pathologies may not occur, but something 
a bit less dramatic may happen which would still make it difficult to reconstruct the higher- 
dimensional cells of the complex from intrinsic properties of the Voronoi graph alone, as we 
will see with a specific example in the next section. 



IV. THE INVERSE PROBLEM 

In this section we discuss how to determine when a. {D + l)-valent graph is a tiling graph 
for a PL-manifold. We first give a procedure for constructing a cell complex from the graph, 
where for greater clarity the steps are introduced for the two-dimensional case, and then 
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FIG. 2: Example with small manifold characteristic scales. The dashed lines are Voronoi edges. 
Some Voronoi "cells" in the handle are not cells in the topological sense, and the edge graph is not 
connected. A similar situation can arise in a topologically trivial, but highly curved M. 

generalized to three or more dimensions. We then prove that those cell complexes are PL- 
manifolds, and discuss the extent to which one in fact obtains differentiable and Riemannian 
manifolds. 



A. Two-Dimensional Cell Complexes 

In two dimensions, the problem is to establish whether a given 3-valent graph 7 is a tiling 
graph, i.e., whether there is a cell complex Q, homeomorphic to a 2-manifold, whose vertices 
and edges are the graph 7. In other words, we need to establish whether there is a way to 
identify sets of edges in 7 that can be "filled in" to become the 2-faces of a cell complex 
satisfying the cell incidence relations. For a set of edges in 7 to be a potential 2-cell boundary, 
it must first of all be closed, so 

(1) Define a (non self-intersecting) loop to be a chain of consecutive edges 

a = {ei, 62, • • ■ , ex}, (4.1) 

where each joins two vertices Vk and ffc+i, that closes on itself and no two vertices along 
it coincide, except for vi = vk+i- 

Most loops are not to be thought of as boundaries of 2-cells. To tentatively pick the ones 
that are, we look for the "small ones" in the following sense: 

(2) A loop a will be called a 2-cell or plaquette if, for every pair of vertices v and v' in a, 
the shortest path in 7 between v and v' is contained in a. 

Once we have found a tentative complete set of 2-cells, we have a candidate cell complex Q. 
This complex is the desired one if its vertices, edges and 2-cells satisfy the Voronoi incidence 



properties (3.1), which will ensure its equivalence to a PL-manifold. 



(3) The cell complex Q consisting of the original graph and the 2-cells defined above is 
homeomorphic to a 2-manifold if every edge is shared by exactly two 2-cells (the remaining 
condition, that every vertex is shared by exactly three 2-cells, then follows). 
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FIG. 3: Failure of manifold construction procedure. The two horizontal lines are identified with 
each other. The rectangles are 2-cells, as are the hexagons. The identification of the octagons 
with 2-cells is ambiguous, however, because there are pairs of vertices in each of them that are 3 
or 4 edges apart both along the octagon and along a loop formed by edges that wrap vertically 
around the cylinder. Instead, the latter loops are identified with 2-cells according to our rules. The 



resulting set of cells does not satisfy the incidence relations (3.1 ) both if we accept the octagons as 
2-cells and if we don't. 

If the cell complex fails the test in step 3, we consider the graph 7 to be non-embeddable in 
a manifold. Notice that 7 may in fact be a tiling graph, but one whose manifold counterpart 
has structure on scales of the order of the cell size or smaller, as mentioned in section III (see 
figure |3|. We take the point of view that even in that case the graph is not interesting for 
a description of semiclassical spacetime, although our procedure may be extendible to one 
that is able to recognize and handle this type of situation in general. One reason we believe 
that this is possible is that the results quoted in section |Tl] on graphs and polytopes seem 
to support the possibility that in a large class of situations there is a unique identification 
of 2-cells with closed paths that leads to the right incidence relations for a cell complex. 
However, until a result to this effect is proved, we will assume that our procedure does not 
lead to the correct construction of a manifold when the latter would have small length scales, 
either from its topology or from its curvature. 



B. Three- and Higher-Dimensional Cell Complexes 

We now extend our construction to more than two spatial dimensions. Again, suppose we 
are given a. {D + l)-valent abstract graph 7; our task is to establish whether it is a Z)-tiling 
graph, i.e., whether there is a cell complex Q, homeomorphic to a D-dimensional manifold, 
whose edge graph is the graph 7. In practice, we need to establish whether there is a way 
to identify sets of edges in 7 that can be "filled in" to become higher-dimensional faces. 

Since the pattern is the same in various dimensionalities, we just state it explicitly for 
D = 3. We begin by finding 2-cells as in D = 2, and verifying their incidence relations: 

(1) Define a (non self-intersecting) loop to be a chain of consecutive edges 

a = {ei, 62, ■ ■ ■ , ex} , (4.2) 
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where each joins two vertices Vk and Vk+i, that closes on itself and no two vertices along 
it coincide, except for vi = vk+i- 



(2) Call 2-cell or plaquette any loop a such that, given any two vertices v and v' in a, the 
shortest path in 7 between v and v' is contained in a. 

(3) The candidate set of 2-cells has the right incidence relations with lower-dimensional cells 
if every edge in 7 is shared by exactly three 2-cells, and every vertex by exactly six 2-cells.^ 

If the 2-cells do not satisfy the requirements in step (3), we consider 7 to be non-embeddable 
in 3 dimensions and there is no need to proceed. If they do, we find a tentative set of 3-cells 
in a similar way, and check the appropriate incidence relations: 

(4) Define a closed set of 2-cells to be a finite collection 

C = {ai, a2, ■ ■ ■ , am} , (4.3) 

where each at is a 2-cell, such that every edge is shared by exactly two 2-cells and aj] this 
makes C homeomorphic to a 2-manifold, as in the 2D construction above. Notice however 
that C may not be homeomorphic to a 2-sphere; whether it is or not can be determined by 

calculating its Euler number x(C) = A^vertices - A^edges + A^plaquettes- 

(5) Call 3-cell any closed set C of 2-cells which satisfies a "cell convexity" property analogous 
to the one in step 2 for two dimensions, i.e., such that, for any two vertices v and v' in C, 
the shortest path in 7 between v and v' is contained in C, and which does not contain the 
set of vertices and edges of a smaller 3-cell. 

(6) A collection of 3-cells defined in this way has the right incidence relations if each plaquette 
or 2-cell is shared by exactly two 3-cells, each edge by exactly three 3-cells, and each vertex 
by four 3-cells. 

For each additional dimension one would add the corresponding set of three steps, up to 
Z)-cells if the graph 7 is (D + l)-valent; the pattern is always the same, and one just needs 
to make sure that the incidence relations one imposes are the ones that follow from equation 



(3.1). Just as in 2 dimensions, a graph 7 is considered to be embeddable in D dimensions if 



the steps above produce cells all of which meet the incidence relations. 



C. Discretized Manifolds 



In this subsection, we show that the cell complexes constructed with the above procedure 
really are topological manifolds; in fact, one can show directly that they are PL-manifolds. 



^ The edge graph of a 2-dimensional square lattice (either infinite, or finite with periodic identifications) is 
4-valent, and one could try using our construction to find out whether it leads to a 3-dimensional manifold. 
One would then sec that each of the usual square plaquettes is a plaquette in our sense as well, but each 
edge is only shared by two plaquettes and each vertex by four of them, so one would not associate any 
3-nianifold to this graph with our procedure. 
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The main tool for this proof is the well-known result [21] that (the polyhedron of) a D- 
dimensional simplicial complex is a PL-manifold if the link of every vertex in the complex 
is topologically a — l)-sphere, where the link is defined as follows: Given a vertex f in a 
simplicial complex, consider the set of all simplices cTj which have v on their boundary; then 
the link of v is the union of all other simplices on the boundary of those cij which do not 
contain v (the simplices are taken to be closed). Essentially, what this is telling us is that it 
is sufficient to prove that each point has a neighborhod in the complex that is homeomorphic 
to a ball. With the constructions of the previous subsections, we have built abstract cell 
complexes, in which each cell is identified just by the vertices on its boundary and the 
structure of the complex by the cell incidence relations. Here, however, it will sometimes be 
convenient to think of the cells as actual topological balls in various dimensions. 

As a first step in applying this result, we need to show how to produce a simplicial complex 
that is PL-equivalent to a cell complex obtained as a result of our construction. Given a 
cell complex f2, there are actually various ways of obtaining such a simplicial complex; the 
one we use can be called a "barycentric decomposition" (although, strictly speaking, we do 
not know which point in the interior of a cell is its barycenter), and is a prescription for 
subdividing all cells in order of dimensionality, starting from the lowest. Edges are already 
1-dimensional simplices and don't need to be subdivided. Each 2-cell is subdivided by adding 
a new vertex (to be thought of as a point in its interior), joining it by edges to all previous 
vertices of the 2-cell, and calling each triple formed by the new point and two adjacent ones 
on the boundary a 2-simplex; finally, the original 2-cell is omitted from the complex and 
replaced by the new set of simplices. In D > 2 dimensions, the prescription then proceeds 
recursively: after having subdivided all fc-cells, each {k + l)-cell, if any, is subdivided by 
adding a new vertex, joining it by edges to all vertices on the boundary of the cell (including 
the ones that were added in previous steps), and calling each set formed by the new point 
and a (fc — l)-simplex on the boundary a fc-simplex; as before, the new set of simplices then 
replaces the original cell in the complex. 

We now need to show that the link of every vertex in this simplicial complex, both the 
original cell complex vertices and the new ones, is a topological [D — l)-sphere. For the last 
set of added points, the proof is simple: the link of each one of them is the boundary of the 
original D-ceW it was added to, which is of course a, {D — l)-sphere. For the previous set of 
added vertices, the "barycenters" of the original {D — l)-cells, one can argue as follows. The 
simplices of the barycentric decomposition that share the vertex in question and are on one 
side of its {D — l)-cell form a cone on that {D — l)-cell; the two cones on opposite sides are 
joined at their common base and their union is topologically a D-ball, with a. {D — l)-sphere 
for boundary. If we limit ourselves to D = 3, the only remaining vertices of the barycentric 
decomposition are those of the original cell complex. For these, the simplest way to show 
that they have a neighborhood in the cell complex that is homeomorphic to a D-ball is to 
use a different construction. Consider one such vertex wq, and the D + 1 vertices (f i, f^+i) 
that are connected to fo by edges. Any set (fi, ...,vd+i) of D vertices chosen among 
the above (the hat means that we have removed fj), together with fo, defines a D-simplex 
Aj. The union U^'^^Aj of these simplices is another D-simplex A, with vertices (f i, vd+i)- 
Then the {Aj} form a simplicial complex in which the link of Vq is the boundary of A; since 
A is a D-simplex, this boundary is a (D — l)-sphere. 
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V. OUTLOOK 



In this paper we have provided a partial solution to the graph version of the inverse problem 
for discrete geometry: When can one obtain a differentiable manifold from a graph? We 
have defined what we mean by obtaining a manifold from a graph, and argued that the most 
natural setting for addressing this question within the context of loop quantum gravity is 
that of Voronoi complexes. We then provided a procedure for constructing, in principle, 
the manifold for a graph of the Voronoi type. The procedure consisted in defining, using 
the graph, cells of increasing dimensionality. Under certain conditions, the resulting cells 
define a cell complex Q that is homeomorphic to a uniquely defined PL-manifold M (up 
to isomorphisms) — the manifold is unique if Q is; but Q is obviously unique, since it was 
constructed using only intrinsic properties of 7, with no arbitrary choices. On the other 
hand, given a generic set of points in a Riemannian manifold {M,gab), distributed with a 
high enough density that the manifold only has characteristic lengths on scales larger than 
the point spacing, the Voronoi construction would give a cell complex on which our procedure 
would work, and would therefore give us back the original manifold. 

To provide a more complete solution to the inverse problem, several limitations need to 
be addressed. The first limitation is related to the types of situations in which the procedure 
does not work. Although, for the reasons discussed in section [ll| we consider (D + l)-valent 
graphs to be the most important ones for semiclassical loop quantum gravity, we do not 
expect them to be the only ones that play a role in that sector of the theory; in fact we 
expect that, to properly take into account quantum fluctuations and non-local features, non- 
manifoldlike discrete structures need to be considered. On the other hand, in practice the 
most important configurations are likely to be the ones that are close to being manifoldlike 
in some sense, and it may be possible to characterize these in a more precise way. From 
the arguments in section [TT| we assume that, if associating a manifold to a graph means 
determining whether the graph is a tiling graph, then generically the corresponding cell 
complex will have the incidence relations of a Voronoi complex; our considerations will not 
apply to other possible ways of finding manifolds from graphs. 

There are two kinds of underlying reasons why a graph may not be manifoldlike in this 
sense: either it does not have a definite valence to begin with, or it does but the higher- 
dimensional cells one constructs as described above do not satisfy the right incidence rela- 
tions. Keeping this in mind, we can identify some types of "almost manifoldlike" graphs. 

The simplest type of obstruction is the one in which the resulting complex can be inter- 
preted as being homeomorphic to a manifold with boundary. In the case of a (D + l)-valent 
graph, this happens when all cells of dimensionality n < D obtained following the steps 
described above satisfy the incidence relations (3.1 ), but there are one or more closed sets of 



{D — l)-cells in which every [D — l)-cell is on the boundary of only one D-ceW. (see figure |4|). 

As a different type of example, consider a case in which the obstruction is "small", i.e., 
there is a finite and small (in terms of what is considered to be macroscopic in the graph) 
number of vertices and edges which give rise to cells that don't satisfy the incidence relations, 
but are such that, if they are removed from 7, the rest of the graph does look like a manifold, 
although this time with a boundary (see figure [s]). Then one has the option of giving up 
on describing that portion of the graph in detail, and replacing the removed set of vertices 
and edges by a single D-cell, to which one can assign an effective volume from the number 
of deleted vertices and the structure of the deleted edges, or of looking at the obstruction in 
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FIG. 4: Graphs and boundaries. According to our definition, the left-hand graph is a tiling graphs 
for the 2-sphere, since even the outer polygon is a 2-cell. The right-hand graph is a tiling graph for 
the 2-disk, a manifold with boundary. 




FIG. 5: Localized obstruction. As it stands, there are too many "2-cells" in the central portion of 
the graph, in terms of their incidence relations, but if the two crossed edges are removed, the rest 
of the graph tiles a manifold with an boundary which can be removed by adding a single 2-cell. 

more detail. 

In addition to genuine obstructions to the existence of an appropriate cell complex for 
a given graph, another open issue is whether one can modify our procedure to address the 
situations in which a {D + l)-valent graph is a tiling graph but the procedure described 
above fails to produce a manifold. 

Two aspects of the subject, the existence and uniqueness of a differentiable structure 
on the manifold we construct and the existence and amount of freedom in the choice of a 
Riemannian metric thereon (the metric is far from being unique), go beyond what we have 
discussed so far. If we restrict ourselves to 2-dimensional and 3-dimensional manifolds, then 
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well-known results from differential topology (see, for example, reference [2T]) guarantee that 
each PL-manifold has a unique differentiable structure on it. 

As mentioned in the introduction, our goal is to explore the possibility of recovering within 
loop quantum gravity the full set of geometrical structures needed to describe a classical 
geometry. In this paper we addressed the construction of the underlying manifold from one 
of the graphs used to define a state, while metric information on the spatial geometry and 
its time derivative will be encoded in the details of the wave function. This wave function 
could be seen either as a state on a three-dimensional hypersurface, or as a 'boundary state' 
for which an amplitude is computed using spin foams ^22j. In this case, the graph has to be 
tetra-valent since it has to match with a cell decomposition on which the spin foam is defined. 
We hope that the construction here outlined will be useful in the quest for a semiclassical 
limit of loop quantum gravity. 



Acknowledgements 



This research was supported in part by Perimeter Institute for Theoretical Physics, by the 
grants CONACyT (Mexico) U47857-F and NSF PHY04-56913, and by the Eberly Research 
Funds of Penn State. 



[1] A Ashtekar and J Lewandowski 2004 Background independent quantum gravity: A status 
report Class. Quant. Grav. 21 R53 (arXiv : gr-qc/0404018 ) ; 
L Smolin 2004 An invitation to loop quantum gravity arXiv : hep-th/0408048j 
T Thiemann 2002 Lectures on loop quantum gravity ar Xiv : gr-qc/0210094| 

[2] C Rovelli 2004 Quantum Gravity (Cambridge University Press) 

[3] T. Thiemann 2007 Modern Ganonical Quantum General Relativity (Cambridge University 
Press) 

[4] C J Isham 1995 Structural issues in quantum gravity 'arXiv : gr-qc/9510063' 

C J Isham 1993 Prima facie questions in quantum gravity arXiv :gr -qc/931 0031 
[5] L Smolin 2005 The case for background independence arXiv:hep-th/0507235l 
[6] K Giesel and T Thiemann 2007 Algebraic quantum gravity (AQG) I. Conceptual setup 

Glass. Quantum Grav. 24 2465-2497 (arXiv : gr-qc/0607099 ) 
[7] See, e.g., A P Gentle 2002 Regge calculus: A unique tool for numerical relativity Gen. Rel. 

Grav. 34 71 (arXiv:gr-qc/0408006) 
[8] See, e.g., R Loll 1998 Discrete approaches to quantum gravity in four dimensions Living Rev. 

Rel. 1 13 darXiv : gr-qc/9805049| ) ; 

J Ambj0rn et al 2006 The universe from scratch Gontemp. Phys. 47 103-117 
( arXiv : hep-th/05090 10 ) 
[9] L Bombelli, J Lee, D Meyer and R D Sorkin 1987 Space-time as a causal set Phys. Rev. Lett. 
59 521-524; R D Sorkin 2003 Causal sets: Discrete gravity (Notes for the Valdivia Summer 
School) |arXiv : gr-qc/0 309009 
[10] For recent reviews, see F Dowker 2005 Causal sets and the deep structure of spacetime 100 
Years of Relativity - Space-time Structure: Einstein and Beyond ed A Ashtekar (World 
Scientific) (arXiv : gr-qc/0508109J ; 



15 



J Henson 2009 The causal set approach to quantum gravity Approaches to Quantum. Gravity: 
Toward a New Understanding of Space, Time and Matter ed D Oriti (Cambridge University 
Press) (arXiv:gr-qc/0601121) 
[11] S Major, D Rideout and S Surya 2009 Stable homology as an indicator of manifoldlikeness in 



causal set theory arXiv:0902 . 0434 



[12] A V Evako 1994 Dimension on discrete spaces Int. J. Theor. Phys. 33 1553-1568 
darXiv : gr-qc/9402035 ) ; 

M Requardt 1998 Cellular networks as models for Planck-scale physics J. Phys. A 31 
7997-8021 (arXiv:hep-t h/980 6135); M Requardt 2006 The continuum limit of discrete 
geometries Int. J. Geom. Meth. Mod. Phys. 3 285-314 ( |arXiv:nia th-phy/0507017|) 

[13] A Ashtekar and J Lewandowski 1997 Quantum theory of geometry. II: Volume operators 
Adv. Theor. Math. Phys. 1 388-429 (arXiv:gr-qc/9711031) 

[14] N H Christ, R Friedberg and T D Lee 1982 Random lattice field theory: General formulation 
Nucl. Phys. B202 89-125; T D Lee 1985 Discrete mechanics How Far Are We from the 
Gauge Forces (XXI Ettore Majorana School, Erice 1983) ed A Zichichi (Plenum) 

[15] E Tonti 1995 On the Geometrical Structure of Electromagnetism Gravitation, 

Electromagnetism and Geometrical Structures (for the 80th birthday of A Lichnerowicz) ed 
G Ferrarese (Bologna: Pitagora Editrice) 

[16] F Markopoulou and L Smolin 2007 Disordered locality in loop quantum gravity states Glass. 
Quantum Grav. 24 3813-3823 (|arXiv : gr-qc/0702044l) 



[17] H Finkel 2006 Stochastic evolution of graphs using local moves arXiv:hep-th/0601163 



[18] L Bombelli, A Corichi and O Winkler 2005 Semiclassical quantum gravity: Statistics of 

combinatorial Riemannian geometries Annalen Phys. 14 499-519 ( |arXiv : gr-qc/ 0409006] ) 
[19] B Griinbaum 2003 Gonvex Polytopes 2nd edition (Springer) 

[20] L Bombelli 2002 Statistical geometry of random weave states Proceedings of the Ninth Marcel 
Grossmann Meeting on General Relativity eds V G Gurzadyan, R T Jantzen and R Ruffini, 
World Scientific (arXiv : gr-qc/0101080 ) 
[21] W P Thurston 1997 Three- Dimensional Geometry and Topology (Princeton University Press) 
[22] See for instance: F Conrady and L Freidel 2008 On the semiclassical limit of 4d spin foam 
models Phys. Rev. D 78 104023 ( [arXiv: 0809 .2280 ) 



16 



